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We propose a method to measure the hadronic matter compressibility by means of the event- 
00 \ by-event analysis of heavy-ion colhsions at high energies. The method, which utilizes the ther- 

0^ ■ modynamical relation between the compressibility and the particle number fluctuations, requires a 

0^ ' simultaneous measurement of the particle source size, temperature and particle multiplicity. 

T-H ■ 
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ff^ ' Large acceptance detectors allows one for a detailed analysis of individual collisions of heavy-ions at high-energies. 

O^ , Due to hundreds or even thousands of particles produced in these collisions, variety of powerful statistical methods can 
be applied. Then, such event-by-event studies can provide valuable dynamical information which is otherwise hardly 
►^ , available. For example, we have shown in |l||, see also [^1, that the correlation between particle multiplicity and their 
f^ ■ average transverse momentum, which is observed in proton-proton interactions, causes sizeable fluctuations of the total 
(T^ , transverse momentum of particles produced in a single nucleus-nucleus collision if such a collision is a superposition of 
^^ ' independent nucleon-nucleon interactions. The fluctuations from Pb-Pb collisions at 158 GeV per nucleon studied in 
^^ , NA49 experiment have appeared to be noticeably smaller than the properly normalized fluctuations from the proton- 
proton interactions at the same collision energy [^ . Thus, a substantial role of the secondary interactions in heavy- ion 
Qv , collisions has been proven without a reference to any collision model. 

Stodolsky Q and Shuryak [g] have made another interesting proposal of the event-by-event analysis. They have 

adopted a standard assumption that the hadron matter from the nuclear collisions is in thermodynamical equilibrium. 

Then, according to the well known thermodynamical relation, the temperature fluctuations can be treated as a measure 

O , of the heat capacity of the hadronic matter. Shuryak ||^ has also briefly considered the relation, which couples the 



I 



:3 



particle number fluctuations to the particle number derivative with respect to the chemical potential. In this paper 
we follow a similar line of reasoning and present a method to determine the hadronic matter compressibility via the 

^ multiplicity fluctuations. 
►v> The method is based on the thermodynamical relation, see e.g. lq,0|, which expresses the dispersion of the particle 

^ number A^ in a volume V through the isothermal compressibility as 



{Nf _ V^ f dp 



{N^) - (7V}2 T \dV 



(1) 

T,{N) 



where T is the temperature and p the pressure. In principle, the relation ffl) allows one to find dp/dV as a function 
of (iV), V and T and then to reconstruct the equation of state. Below we discuss how to realize this program in 
heavy-ion reactions. 

There are many sorts of hadrons (7r+, tt", tt^, K'^, K^ , ) in the final state of nuclear collisions. Therefore, 

our first task is to generalize eq. ffl) to the system of k components. To solve this problem one has to assume that 
the hadron system at freeze-out is not only in the thermal but in the chemical equilibrium as well. The distinction 
between the thermal and the chemical freeze-out is discussed at the end of this paper. The extensive analysis ^^ of 
the experimental data from AGS and SPS accelerators shows that the assumption of the chemical equilibrium is very 
well satisfied in heavy-ion collisions at least for nonstrange particles. 
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To derive the relation analogous to (||) we follow |^ and write the grand canonical sum for the fc— component system 
as 

E{zi,...Zk,V,T) = ^ z^^.-zj^" Qni,...nAViT) , 

Nu...Nk 

where Zi is the fugacity of the i— th component and Qn^,... Nki^^T) is the canonical partition function. Keeping in 
mind that 

with (3 = T~^ and F{Ni,... Nk,V,T) being the system free energy, one finds that the probability PNi,...Nk to find 
{A^i, ... Nk} particles is proportional to 

In the case of a hadronic system the particle numbers {Ni, ... Nk} are not conserved. Nevertheless the free energy 
F{Ni, ... Nk, V, T) with the numbers of particles {iVi, ... Nk} being fixed is still of physical meaning. However, this is 
F{{Ni), ... {Nk),V,T), which corresponds to the thermodynamical (thermal and chemical) equilibrium. The average 
particle numbers {(Ni), ... (Nk)} are found as a minimum of F{Ni, ... Nk, V, T). If the particles carry charges, say + 
and — , and the total charge Q is conserved, then the average values {N^) and {N~) are also found as a minimum of 
F but the additional constraint {N^) — {N^) — Q must be imposed. 

When we deal with the particles obeying quantum statistics, the approach described above, where the average 
particle numbers arc found as a minimum of F{Ni, ... Nk, V,T) is rather inconvenient because there is no compact 
expression for F{Ni, ... Nk, V,T) even for the ideal gas. Then, we introduce the chemical potentials. However, this is 
a technical problem and not the matter of principles. 

After this comment we expand the free energy around the average particle numbers and take into account only the 
first three terms i.e. 



FiN,,...Nk,V,T) = F{{N,),...{Nk),V,T) 



(2) 
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dF{Ni,...Nk,V,T) 



dNi 



iN^ ~ {N,)) 



Ni = {Ni).,...N„ = {Nk) 



l^&'F{Ni,...Nk,V,T) 
9 ^ 
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dN.dNj 



(iV,-(iV,))(7V,-(7V,)) 



Nt = (Ni),...Nk = (Nk) 



Due to the relations 
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A,idF{Ni,...Nk,V,T) 
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where fii is the chemical potential of i— the component, one finds 



Zi = e 



/3m> 



Pnu-.-N^ --exp 



Ij2a,,{n,-{nmn,-{n,)) 



«j 



(3) 



with the matrix A given as 



Ay = (3 



d^FiNi,...Nk,V,T) 



dN.dNj 



As well known, see e.g. pO|], the moment matrix Mij of the normal distribution (m equals 

M,, ^ ((TV, - {N,)){N, - {N,))) = {A-% . 



(4) 



(5) 



We are now going to relate the matrix A to the compressibihty. Since the free energy is an extensive quantity, 
which can be expressed as 

F{N,,...Nk,V,T)^Vf{pi,...puT), 

with Pi = Ni/V, one shows that 

dcf dF{Ni,...Nk,V,T) 



piNi,...Nk,V,T) 



dV 



FiNi,...Nk,V,T) , ^ dF{Ni,...Nk,V,T) 
V ^l^P- dN, 



and then 



dpiNi,...Nk,V,T) ^ d''F{Nu-Nk,V,T) 



dV A^^'^^ dN.dN, 



Using eqs. (ph and (P) we get the final result 

~^^- =- 

T \dV 



J2{N,){N,)iM~%=-^(^) . (7) 



hJ 



T,{N)i,...{N)k, 



The relation M) essentially simplifies when the system components are independent from each other. Then, the 
matrices A (H) and M (|f) are diagonal (My = A^ = for i ^ j) and 

V im' V' fdp\ 

The question whether {{Ni — {Ni)){Nj — (Nj))) = ior i ^ j can be answered experimentally. 

As already mentioned it is somewhat unusual in the statistical hadron physics to consider the thermodynamic 
quantities, like the compressibility from eqs. (|7|i|8|), at the fixed hadron average numbers. Computation of such 
quantities is straightforward if the hadrons are assumed to obey the Boltzmann statistics. In fact, this is a quite good 
approximation at the freeze-out stage when the hadron gas is expected to be rather dilute. 

We consider as an example the classical multi-component van der Waals gas. The equation of state is taken in the 
form 



p(y-^(iV,:K:)=^(Ar,)T, 



J 



where the parameter Vi is related the volume of the i— sort hadron. The r.h.s. of eq. (R) or (@) then equals 



TAN)r....{N)„ (^l_J2,m^^ V ^ 



where the last approximate equality assumes smallness of the van der Waals correction. When u^ ^ we get the 
ideal gas limit. The compressibility of the classical ideal gas is obtained from the relation ffl) if the moment matrix 
M (I) equals Mij = 6'^N,). 

When the particles obey quantum statistics one usually introduces the chemical potentials to compute the ther- 

modynamical quantities. Then, the pressure is a function of {pi, pk} (not of {(-/Vi), ... (Nk)}) and we express the 

derivative at fixed {(A^i), ... {Nk}}, which is present in eqs. (0,||), through the derivative at fixed {pi, pk}- This is 

done by means of the respective jacobians, see e.g. |0]. We give here the formula for a single component system 



dp\ _ ( dp\ \9t^)T.y\^^ /T,M 



dyJr.XN) V5^/T,M 



fdp\ (aW\ 



(d{N)\ 



and apply it to the ideal gas of pions which are bosons with three internal degrees of freedom. Then, 

p(M, 1/, T) = -ST / ^ ln(l - e-^f^-'^)) , 



where Ep = ^vn? + p^ with ra being the pion mass and p its momentum. For massless pions with /i = 0, which 
corresponds to the chemical equilibrium, we get the r.h.s of eq. (pi) in the form 



V^ {dv\ _ 18(2(3) ^^^3 _ 6C(3) 

T \dY 



VT^ = ::i^ {N) ^ 0.71 {N) , (10) 

T,{N) '^' 



where C,{z) is the Riemann function and C(3) = 1.202. Comparing eqs. (^ and (10) one sees that the particle number 
fluctuations are larger in the ideal gas of bosons than in the classical one. 

Let us now discuss how to use eq. (0) or (g) to determine the hadron matter compressibility. The l.h.s. of these 
equations is fully determined by the particle multiplicity but to extract dp/dV one also needs to deduce T and V in 
individual collisions. It is not a trivial task to get these three quantities. 

One should first realize that the final state characteristics correspond to the so-called thermal freeze-out when the 
whole system disintegrates due to the switching-off the inter-particle interactions. However, the system chemical 
composition is fixed earlier, at the so-called chemical freeze-out when the interactions, which change the particle 
numbers, are no longer effective. The particle multiplicity, the system volume and temperature, which enter eqs. MM) 
refer to the chemical freeze-out. Therefore, these quantities, which are obtained from the final state characteristics and 
consequently correspond to the thermal freeze-out, should be further recalculated to get their values at the chemical 
freeze-out. One can refer here to the procedure described in e.g. P], where the system volume and temperature at 
the chemical freeze-out are inferred from the multiplicities of the different sort particles. 

The second but related difficulty lies in the fact that a sizeable fraction of the final state particles originate from 
the decays of hadron resonances. Since the particle multiplicity from eqs. mM) is that at the chemical freeze-out, the 
final state multiplicity should be recalculated to reconstruct the resonance contribution. This can be achieved within 
the thermodynamical model which takes into account hadron resonances, see e.g. MM- 

The system temperature at the thermal freeze-out is usually found as an inverse slope of the transverse mass 
spectrum {rriT = \Jvn?- -|- p^ with pT being the particle transverse momentum). This is however the 'effective' 
temperature which incorporates the effect of transverse collective motion of the hadronic matter pT| and that of 
the resonance decays |12| . The actual temperature can be disentangled by the simultaneous analysis of different sort 
particles ||ll|.|l2[. The procedure is not very exact but the temperature is presumably measurable within 10% accuracy. 

The system volume in nuclear reactions is controlled due to the collision centrality selection. In this way one can 
change the volume in a rather broad range. One usually deduces the particle source size at the thermal freeze-out by 
means of the particle interferometry measurements, see e.g. |lj]. The high multiplicity of relativistic nuclear collisions 
allows one for the event-by-event interferometry analysis. However, one should keep in mind that the hadron system 
decay is a dynamic process and the concept of the system volume is then not very precise. The hadron resonances 
also complicate the analysis, see e.g. ||lj]. In spite of these difficulties, we still believe that the system volume can be 
measured with an accuracy which enables one to extract the hadron matter compressibility from eq. (M). 

How the data should be collected and processed? First of all one should eliminate or at least reduce the trivial 
fiuctuations due to the impact parameter variation. This can be achieved by the trigger condition. The collision 
centrality is well known to be strongly correlated to the transverse or forward energy observed in the collision. 
Therefore, selecting the transverse or/and forward energy from a narrow interval we can collect a sample of collisions 
of a similar geometry. Then, one determines the system volume, temperature and hadron multiplicity for every event. 
The events are split into subsamples in such a way that a subsample contains the events with coinciding or close 
temperatures and volumes. For every subsample one constructs the multiplicity distribution and computes the l.h.s. 
of eq. (n) or (H). Then, one gets the compressibility for a given V and T. Combining the results obtained for several 
subsamples we get the compressibility as a function of volume and temperature. The range of Y and T will be 
presumably rather small. To get dp/dV in a broader domain one should vary the collision geometry by means of the 
trigger, change the projectile-target system and the collision energy. 

In summary, we have proposed a method to determine the hadronic matter compressibility due to the event-by- 
event analysis of the particle multiplicity fluctuations in heavy-ion collisions. The thermodynamical relation, which 



connects the two quantities in the multi-component system, has been derived. As an illustration we have considered 
the compressibility of the classical van der Waals gas and the quantum gas of massless pions. Since the method 
requires not only the measurement of the particle multiplicity but of the system volume and its temperature as well, 
the problems related to such measurements have been briefly discussed. Finally, the procedure of data analysis has 
been suggested. 

I am very grateful to Marek Gazdzicki, Mark I. Gorcnstein and Edward V. Shuryak for critical reading of the 
manuscript and fruitful discussions. 
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